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ABSTRACT 


of flow 

; deviations fron Darcy ty^''/-ic obcorverl in many 
field siTno L 3 011.0 . ao such it has become necessary xo 
use liecd loos ecjuationo su^'gesxed by Forchheimer cr< 
rio'boch for realnsbic analysis of such problems., 'A 
review/ of bhe literature on the subject of ITon-linear 
Plow Through J orouo Media is presentorl. Analytical 
sulubionc to some of the non-Darcy type of flow problems 
ore 2>resontecl based on Du pits assumption. The effect 
of nonlinearity parameter on discharge ratio and draw- 
down GurTG are discussed in detail. Equation.?! are 
* 

developed for nonlinear flow through nonhomogenous 
porous Media, ^ 

^'She experimental inves i^igotions of nonlinear 
coeiXiCient.c a and b for various materials are presented 
A method of i^redicting the values of a and b using the 
ovs'iJoble data is suggested. Equations developed for 
nonlinear flow through nonhomogenous porous media are 
verafied experimentally for a two layer gravel bed. The 
estimation of head lo.ss through gravel pack for a tube 
v/ell ir. illustrated, 



CIlArT^R 


I 


Tl'TORCDUCIirr 


1 . i I'he problems involvuifi flo\/ of '/i^ater through porous 
]i' 0 ( \c vovy v/iaely ixi scope end its engineering applications, 
frefli licnally these problems have boon solved based on 
harc^^'s linear relationship between head loss and 
vflociu/ ; 

V = Ki (1 .1) 


\/herc l~ = Darcy's coefficient of permeability, 

dH. 

1 = = negative hydraulic gradient 

dS 

Hg, = Total heed i,e. sum of the piezometric 
and velocity heads, 

S = Distance measured in the direction of 
ve] ocity 

Y = Bull: or average seepage velocity. 


The equation (1.1) satisfactorily describes the 
flow conditions, provided velocities are small i,e, the 
flow 16 laminar. However, as the velocity of flow, 


pa.rtxcle size and Reynolds number increase Darcy's linear 
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rclabiopship leads bo iiiaccurotp prefiictions . Thus while 
ii’atiy piJctical problems of floiv through por^'iis materials 
r U' he solver’ accuiatelv based on Darcy's lav/, various 
iielr" situations have been ol^scrveol v/here a more accurate 
ro] a bionship bet’/een heed loss snd velocity must be 
ej.ipj.oyc'd CO obtain realistic .solutions. 

Goife ol the situations where in the use of 
nonliiior'i relationship becomes necessary are ; 

(1) Dlov/ in Che area adjacent to a pumping 
iiell e .peciaily m a coarse grained aquifer. 

(2) Dlo\/ in the filter pechs (gravel packs) of 
lubi" \;ell& . 

(3) PloY/ tlirougn rod' fill banks and dams, 

(ri) Plov/ rlirough filters used in water purification 

pi ontfa , 

In the recent years, many investigators have 
rti^'died ground water flo'" problems involving nonlinear 
ilo’' sjtuaticns. Hov/ever, the analysis i& far from 
coT'iplei.e, as such an attempt has been ruade in this 
investigation to analyse some of the nonlinear floY/ 
problem.^ and to verify the equations developed for 
nonlineor Ilo'/r through nonhomogen^us medium. 
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! ,2 Literature Review 
1.2^1 Re i.-3 1 a ear Rqu at i ons 

Tne nOia-Darcy regiif'e of flov' has beer generally 
•lescribei''' by the enuation suggested bj' Porehlieimer , 

1 ^ aV + b''~ (1.2) 

in 'hiich 0 and b are constants determined by the properties 
of bhe fluid and porous loeuia. Pumber of authors (V’olker~ 
1963), (Ahi.'Cd S Suneda - 1969), (O’Heill, Parkin, I’odd, 
Ranganadna Rao S Suresh - 1970) have supported 
bhis relotionship after conducting exTierimental investigob- 
ions. By deduction from Ravier Gtokes equation ihmed 
and -junada (1969) shov/ed that the governing head loss 
eque bion is of the form suggested by Porehheimer, i.e. 
equation (1.2)., but values of coefficients a and b are 
nob sbricbly constant and depeiid on the Reynolds number 
of bhe flov/. However, the available data indicates that 
a and b are essentially, const tints over a range of Reynolds 
nunbers , 

Another form of equation suggested bjr Missba.ch : 


1 = 


(1.5) 
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IS also in coiiiron us«, . Ixi this equation, C is a 
coefiicieiit de'erijiinoa one pa^opeities of the riuid 
end j.edium. 

m = an exponent lyin^ beL 'een 1.0 and 2.0 The value 

01 iO. vaiies froii 1,0 for Iniuina'' flo’' to 2.0 for complebely 

ourlujent fllo.v. 

^he jse of eciuation (1.3> lo supported by nany 
aubhers. (Anandakrishnan, Vaisdori. .^ala ~ 1963) > (Parkin, 
Trollope, Law-on - 1966). 

Vilkins (1955) j using wide range of aggregate 
sizes and garbles, derived the eoiualion of the form. 

(1.4) 

where Vy = Seepage velocity 
1 

= in 

r^ = Hydraulic mean radius 

_ Void ratio 

Surface area per unix volume 

= Dynanac viscosity of the fluid 


are constants 



5 


The value of , varied from 32.9 for crushed 
stone bo 46.5 for marbles in inch units. 

1.2.2 The Beginning of Fon-linear Blov/ 

The acce’ited jjiethnd of expressing tne limitations 
bo Boxc.'-'b law has, been by pipe anoiogy to determine 
the Ileyni'lds number at the point 'vhore departure from 
tne linear lelationship as first noted. It is v;ell known 
that chc Keynolcis number 

, f!L (1.5) 

\'hcrc p = Tensity of the fluid 

ju = Dynamic viscocity of the fluid 
1 = Characteristic length of flow i .e 

median sa ze of the particle. 

’'hilo calculating Reynclds number some 

inveu Mfot ore (Porkan^ Lavfsonj Tr^^llo c - 1966) used 
seepage velocity v/hich rs function of poiosity (n) . 

Seepage velocity or void velocity (Vy) is 
defined in terms of void ratio 'e'. 



Q 

I 


e 


(1 . 6 ) 



G 


_ Total disch£-r£'e 

” Area of voids 

whenre, A = Total a^ea 

Limiting values of Reynolds number (Rg) from 
0.1 to 75.0 have been measared by various investigators 
(Anandaicrishnan, Varadaiajlu - 1963 ), (’’/right - 1968 ), 

(VolLer - 1969 ) (Ahmed, Gunada ^ 1969 ) for wide range of 
fluids and materials The large variation in the 
Reynolds number (Eg) reflects in the deficiency of the 
usual form of parameter in accounting for the variables 
involved in the flow, such as particle size, shape, 
particle arrangement, roughness and porosity, 

Eventhough the laminar flow tlirough a porous 
media is considered analogous to the laminar flow through 
the p3pe, the analogy does not hold good beyond linear 
regime of flow. The logarichmic plot of resistance 
coefiicient (P) versus Reynolds number (Rg) does not 
show a sharp jump in the resistance coefficient found in the 
the case of pipes at the end of laminar regime . The 
plot, hov/ever exhibits a long progression from the 
stream line Darcy flov. state in which head loss is 
proportional to the velocity, to a turbulent state in which 
the head . loss vanes as square of the velocity (Wright - 1 968 ) , 
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(-‘■ix'.ed, dunada >- 1969) The gradual tra/isition found to 
e::i&fc beb’-een the regions of lemmar and turbulent flow 
i& irooobly due to the large variation in the particle 
&i 2 ie ond voids v/ithin a given porous medium. 

The limitations of pipe flov/ analogy is that in 
^.ipe flow the motion is along rectilinear paths at constant 
velocity for steady state conditioii. VhiJe in case of flow 
through poious media, the fluid undergoes continual cycles 
ol c'ccelerotion and deceleration v/nile follov/ing curvil- 
ine<3r paths ♦ 

The major difficulty in the definition of 
rieynclds number (R ) is the length dimension, for 
pipe flow this dimension is clearly defined by 
pipe diameter, but for porous medium the flow channel 
is so variable in dimensions that a similar dimension 
IS almst 12 /ipossible . The effective particle size 
appears iio be more significant than a flow path dimension# 

Taylor (1948) calculated that for a unit 
hydraulic gradient, the mu,ximum diameter of the uniform 
sand in \7hi ch the laminar llov' wi31 occur is 0.5 m.m. 

On the basis of the experimentol results and 
interpretations of flow through porous media, bright 
(1968) has classified the regimes of flow into four 
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cate^oriea as follows, 

(1) A laminar regime in '"hich at every yjoint 
micrnvelocifcy is &tai:ionery and head less is directly 
pj oportional to the velocity. The viscns forces 
predominate and maximum velociby occurs near the center 
of each flow passage . 

C 

(2) A steady in^erial regime in which at every 

A 

point microvelocity is still stationery but head loss 
hi-ts ceased to vary linearly with velocity. Both 
VISCOUS forces and intenal actions influence the motion. 
Staxionery vortices may be formed at the upper end of the 
regime . 

(3) A turbulent transition regime in which the 
microvelocity fluctuates at any point with increasing 
but regular frequency and head loss approaches dependence 
on the square of the velocity, 

(4) A fully turbulent regime in which all 
parts of the flov/ are turbulent, the microvelocity 
fluctuates randomly about a mean. The head loss is 
close to dependence on the square of the velocity 

but if a viscous sub-layer is present it will continue 
to decrease slightly as the flow increases. 
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Bccd-use of the Stoch^scic nature of flow 

0 f 

0 

chrougli porous media, it is not to be expected that 
these regimes con be defined precisely. Howeve'” 
foliov/mg range of Reynolds number (R^) can br generally 
accepted (’./right - 196S). 

Regime Reynolds number (R^) 

Iiaminar 1.0 - 5»0 

Steady In\e 2 ^ai 90 - 120 

Turbulent Transition 
Transition about 800 

Nonlinear e^iuations can thus be considered 
to be yelld for a region of Reynolds number (H^) 10 - 
2000 in which steady Interial flow occurs, 

1,2.3 Rriction Factor in Porous Media Plow 

Ahmed and Sunada (1969) showed that in porctis 
media flow a unique relationship exists between two 
diraenaioniess parameters which are defined in terms Of 
the physical properties of the fluid » porous medium 
and flow phenomena. 

The equation derived by Ahmed and Sunada 4^ 
of the form, 



1 


(1,7) 
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and 


R 


+ 1 <0 


'( 1 . 8 ) 


J-" 

v/hore = P’pessude £,radient 


<4.nd k = 


O^d' 


(1 .9) 


k = Soecific perii‘30bilit,y 


C 


n 

c 


prof 


d 

G 


= Coefficient depending on the physical 
of the media* 

= Particle size 

= Acceleration due to gravity 

Priction factor (P-j) = gC^d ^ 


The constants a and b are 


respoc Lively . 


fgk 


- -and 


1 


S 




The results of 18 tesLs by Ai-mied and Sunada 
( 1969 ) showed that a single relationship exists between 
friCbion factor (P^) and Reynolds number (Sg) • Ronlinecritj' 
becomes predominant at Rg = 0.2 and when R^ IB large 
H Ojj.n-oaclieb unity. Ho'^^ever, they have not clearly 
rcomoned the significance oi the parameter la the 
eciuation (1 .?). If fche range of the parameter Q >2 
known from the experimental values for diffesnt materials, 
11 would have been useful bo determine the specific 



11 


permeability sccurateiy Siil to check the relationship 
^ivcn by equate on (1.8)^ Todd (195S) hes defined 
friction factor as 

p = (1.10) 

2 o 

2pLV- 



(1 J1) 


= Precsure drop over a length 1. It has 
beexi shovm that tranaition fror laminar to turbulent 
flow occurs when is in the range betv/een 1 and 10, 
chus indicating an upper limit for the validity of Darcy's 
law . 


ivard(l970) has defined friction factor as 




+ 0.55 C, 


1/2 


(1 . 12 ) 


The equation (1.12) is siiular to the equation 

1 /P 

(1,8) except for the term 0.55G^^ where is variable 
depending on the physical properties of the media, 

Ranganadha Rao and Suresh (1970) concluded 
from their experimental investigations that the relation- 
shii given by the equation (l,8) proves only the validity 
of Torchheimers equation (1 ,3) but cannot prove the 
correctness of the values a and b calculated from 
hyd^^aulic measurements. 
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Since friction factor (?) 15 vaiioufely defined by 
many investigators, iXidgeon (1''164), 'fodd (1970), Ahmed 
and Oiinsua (1969), Ward (1970), there is a need for standar- 
dising a form of the friction factor. The friction factor 
(P^ ) veivu,- Reynolds nuiiber fP^) plot obtained by Aliraed 
and [.iinodo (1969) for different ■^latenala exhibits a 
sin^jle Dine rela fcionship. ’’^here a,s Itidgeons plot shows 
different lines for various materials The plot obtained 
by Todd and Tyag:i (1970) ais'"’ does not agree with Ahmed 
and bunadas plot For different porous materials at a 
pa:’tiGuDar Reynolds number (Rg), friction factor (P) 
cannoc be same since it is function of particle size, 
shape parricle arrangement, roughness and porosity. Thus 
iL oppoors the relationshiij bet'^/een fncbion factor (P) 
and ReynoDo's number (Rg) can be oetter represented by 
diffeenb lineF^ for various mateiials, 

1 .2.4 Fon-lineor Coofficienbs - Values of a and b, c and m 

Many inveotigat ors have obtained the values of 
8 and b for differenb mateiials ovei” wide range of Reynolds 
number . The values obtainedi by Dudgeon (1964), Sunada 
ahrusd ( 1969)5 Ranganadha fao and Suresh (1970) are shown 
in fcablei (l), (2) and (3) respectively. Anandakrishnan 

■j 

aiyt Varadara^lu (1963) obtained the values of K' = q and m 
for different grades of sand (Table 4) . 
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Volker U969) cbta: nerl the values of a, and b 
in Forelineij'iers cinauxon end C c-r-i ro in errjjonential 
equation by loborotox''y ’die" cf j'odel gravel banks 
in Oji o^en flume ♦ 

blixl e tabulating mexi v&l'ier for coefficients 
a ond L in the equation (i.T) and G nncl m in tho equation 
(l,3)j many investigators have not irentioned the range 
of Rejmrldo number over v/hich the coeflicients for various 
mctorials are valid. Since the coeflicients are function 
ox ReynrDciO number (Bg)^ porosity and other properties 
of tne media and water, these results are of limited 
use in practical application. Proper sji^ifting of the 
repoited data and its classification is needed to make 
them usefu] in practical applicaci on. One of the aim 
01 the present study is to present a method of predicting 
the vulue of a and b for a given porous media. 

1.2.5 Analytical and Experimental Invest] gat ions 

^'lany investigators (Parkin, Lawson, Trollope-1966) 
(Robin Curtis - 1967), have studied the flow 
thi-ough rock fill banks and analj'sed the stability of 
slope against seepage. Parkin (1963) made extensive 
studies OX- model dams in a 4feet wide flume and foimd 
the method of gradually varied flow is satisfactory to 
oetormnne the water surface profiles. The flow of air 
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tlirouj’i jo-'QUf i-'ipuia .v:.!? r.biJdieJ bj' /'right (1963)"^ Both 
velooxty ond burljulencc r, e^?ui eiiioat& 'ere n-acle with bhe 
use 01 a hoc v/ire aiiencLiOucr, hri3ritf results show that 
althc'iv^^ linear resi et-mc'- lelatiOii ceose to be "''alid 
ac Royiioido munber (Hg) oi about 2.0, velocity fluctuations 
do no'c begin until Feynoicu nu^ribor (Fp) is aboub 100 and 
Lurbolericc in not fully es-toblished i-ntil it lo about 800, 
In another set, the flow of waber through grovel bed in 
cyli.ioricail and converging paranetures was studied and 
resulio 'ere corrpareij . Prom the observations it is 
found fchab the resistance begins lo deviate from Darcy’s 
linear low v/hen Rg = 2.0 and where R^ = 100, resistance 
increases to 3.5 times the laminar value. These obser- 
vations confirm that Darcy's law cease to be valid 
before the onset of t\Tbulence . 

Volker (1969) carried out the expcriiT'ents on 
flo\/ of waoer through nodel gravel banks in an open 
flune. 'fhe experimental results were checked by finite 
elemcnb nebhod. Dudgeon (1967) studied the wall effect 
in permearaeter using cylindrical and box type permeameters^ 
Ldie err cos due to wall effect ranging from 5 to 15^ were 
Pleasured for particle to permeameter tube diameter ratios 
from 1 ; 5 to 1 ; 250, 
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Ealcic (1566,1969) i<^^c extent', iv^ 'it'jdy of 
iiOiilixieOr Ilou problor’c . Ho h. j iio ciiiiecl on electrical 
t-noloivc iiodel fox' tiie &ol'-Lir,.i ox' &uch probltsnis , 

KadLav aajd 6 ubraniaaya ( 1970 ) hcrj-p nrexeaxeu 
ooluGxcn,*- to the &one of the coniici eeepe^e problems by 
usins P'"' '’ehheimer.sand aloo Pjpits osoun'ption. 
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GtliiPTER II 

i-iITALYTICAL SOIUTIOPC TO COIE TLO TR 01 LEMS 

2.01 In this chaptei sGiutionc to oon'o siniple problems 
with non-Darcy flow situations "re analjfsed by e' 4 uations 
().2) and (1.3) and with the assumption, Dupits theory 
to be valid, 

2.1 riov; into a trench - Artesian "Plow 

Digure (1), shows a constrijction of a fully 
penetrating trench in a stratum under an artesian condition 
\/ith all the relevant dimensions. In the development of 
this equation and subsequent equations, it has been assumed 
that Dupit - Dorehheimer assumption (leonards) holds good, 

2.1 .1 Case (1) 

If Q IS the discharge into the trench per unit time, 
rhen average velocity V is, 

“ HD 

where B and D are width and depth of the previous stratum 
respectively. 
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Substituting the value of V xn equation (1,2), 
dh 

1 = 

<iy 


aQ 

BD 


+ b ( 


BD 




( 2 . 1 ) 


\/here U is the height of v<ater above the hot com of 
trench and y distance from the trench. 


Integrating equation (2.1) for h v/ith boundary 
conditions, at y = 0, h - h. 

C? 


h + 

e 



+ b 



y 


( 2 . 2 ) 


Por boundary conditions at y = L, h = H in the equation 
(2,2) and solving for Q 


/ — ^ Q - ( ■) = 0 (2.3) 

L b2J 


where and H are v/ater level on the face of trench and 
at distance L from the trench respectively. 

v/ith Q = ^ and ) 

St 

where is the flov; into the trench for linear case i.e. 
when b - 0, Qji = (H-h^) and is a parameter, 

signifying the nonlinearity parameter b, the equation (2.3} 
can be written into the following dimensionless form. 

H- g; - 1 = 0 

1 


(2.4) 
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in '/hich in the rauio of nonlinear discharge to 

1 ijiear discharge. 

oolutiGu of ol'f eirati on (?.'^) is shc’vn ii figure (4) 
iOL function of '7^ on semlog plot. For increasing 

valnos of ^ i.e. increasing nonlinearity the actual 

Li 

uioChBrge decreases logorithrncallj' compared to the 

discharge estiiuted on the basis of linear Ian, Thus for 

any poious media defined by (a, b) increasing head (H-h^ ) 

do noc proportionately increase discharge and value of 

discharge ratio decrease with increasing heads. Thus 

the effect of using Darcy lav/ (linear lav/) for non-linear 

regime is to overestimate the expected discharge. The 

expression for drawdown curve obtained from equation (2*2) 

after substituting for Q is same as one obtained by 

(H-hg) (I-Y) 

using linear law i,e, H-h = ? 


2.1.2 Profile of the Phreatic Surface 


The equation (2.2) can be written in the following 


f orm . 




(2.5) 


Por unit width of trench, B = 1, the equation (2.5) can 
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bo rep: 0 sent eo in tne follov/in^ form. 






_JL. 


( 2 . 6 ) 


v;here e'luol to g , discharge (diriensionless) 


for nonlinear case and is the nonlinearity parameter 

o 2 

1 .e . b/a"^ . 

q = Discharge per unit v/idth of trench. 


Dor linear case, tne equation (2.6) reduces to 




1 .0 


(2.7) 


where is discharge (di:aensionle&s) for linear case and 

H 

equal to — rr-p 

Ip = Discharge per unit width of trench for linear 
ca se . 

The solution of the equation (2.6) and (2.7) is 
represented in Figure (^) on sepu logarithmic paper for 

V h 

^ versus with varying values of and same value 

of discharge (dimensionless). The effect of nonlinearity 
for same discharge is to increase the drav^down. For 
values of ^ ^ 0.010, the increasing value of has 

significant effect on draw-down curve . 
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2,1,3 Case (ii) ; Solutiorj by being Pov/er lav/ 

In this ca&e solution to the above mentioned 
proolem is obtained by enoation (1.3) 


1 = 


dh 

dy 


Substicuting the valne of V in eqiiation( 13) , 


cyi 


BD ^ 


(2*8) 


Integrating equation (2.3) \/ith boundary conditions 
y = 0, h = hg and y = 1, h = H, 

H-h. 


c ( ^ - ( 


) = 0 


(2.9) 


with oC - Q/Qjj and 


wliorc 0_ =r KBD 


(H-h )' 


H-h^ 

Tc“ 


when m = 1 . 


The equation (2.9) can be written in the follov/ing 
Oimer sionless lorm. 

m iti 


"V — 1 =■ 0 

t 3 


where 'yj is the nonlinearity parameter, 
3 


(2.10) 
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oolutior-s of the eiiuation (2.10), are presented 
in (6) f or •'•cr''- 1 IS 'yj cj ^ernilogarithnic plot. 

s 

'ith lucressing value of i.e. ircreasing nonlineari cy , 

tile actual discharge cl'=>crei&es logariXiirnically , The actual 

discharge increases \'ith increasing value oi' in for values 

of ^ upto 1.0, when m = 1 , linear case, the discharge 

3 

ratio remains same for all values of . i?or value of 
= 1 .0 the effect of nonlinearity on discharge ratio 

lb nil, For low values of *<1 1.0) the actual 

discharge is higher than linear case discharge and vice ve 
versa , 


2,2 Flow into a Trench - Gravity Flo\/ 

2.2.1 Case (l) 

Figure (2) illustrates the problem under 
consideration. 

The average velocity of flow 'V' is given by 


V = 


Q 

B.h 


( 2 . 11 ) 


where B.h is equal to area of flow. 


1 = 


to 


Substituting the equation (2,11) In equation (1^2), 

_ _ ^ , bQ^ 

Bh j2 j,2 


( 2 , 12 ) 
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Inte£,iating equation (2,12) v/ith boundary conditions 
y =0, h = hg, and y = L, h ^ H, 

(h2-h|) - A (h-h^) a2 10 


aQ 

^ y 


( 2 . 13 ) 


where A = 


bQ 

aB 






-oL 


( 2 . 14 ) 


h. 


where drav/down ratio is equal to — § and'^ i& equal 


h H-h^ 
to ( ^) 


(1“Z^) _ 




(1-Z) 

a-i-oi) 


2 2 0+r^i) 




SO-r^i) 


log. 


2 + ) 
2.Z 

1 


(1- y/i) 


(2.15) 


where Z is equal to -g 

Equations (2.14) and (2.15) can be solved by the method 
of iteration for calculating values of and Z respectively. 



Th? solution to the equation (?. 14 ) represented 

in the iigure (7). It can be ceen trorn the figure (7), 

thac bhe vtlue of I etaeen .01 and 1 .0 discharge 

raxio decree- 5e& logoi ithriicnlly end discharge ratio 

inci eases again f or 1 .0 , I'uus the effect of 

i T ^ 

nonlinearity paraireter ) ma^r decrease or increase 

'I 

dischaige ratio depending on the vclue of For seme 

drn\/do\/_i discharge in alv'ays leas for nonlinear flow than 
that lor linear flow. The effect of drawdo\/n. ratio 
on discharge ratio is insignificant and for values 
of \ .10, it IS practically ml, 

2.2,2 Variation of Discharge with Drav/do\/n fiatio 

The equation (2,15) may also be written in the 
folio''7ing form. 


For B = 1 , 



v/here Q ^2 dimensionless discharge and equals to 


( 2 . 
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For ^ ratio equals to 20.0, solutD on,: to the equation 

(2.16) aio shc'n in figure (8) on semi logarithmic paper, 

H 

'Phe effect Ox nonlinearity parameter i.e.^^ on ratio 

2 e 

for a given value of ^ almost same for all values of 

i\) H 

'12 bet’./een .01 to .1u, and for lo\/er value of ^ > h — 

2 2 e 

ratio increases more rapidly than for higher values. 


2.2.3 Case (ii) . Solution by Pov/er Law 

Solution to the unconfined trench flow by using 
equation (1 .3) 

1 = CV^ 


dh 

•iy 


Q \m 




( 2 . 17 ) 


Integrating equation (2.17) v/ith boundary conditions 
0 t y = o, h = hg and y = L, h = H. 

1 ^ mi 1 

irfr— < 2 - 18 ) 

when m 1 , Q = linear cage, 
with ^ 
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1 1. Tri+ 1 

^ ~ ‘‘e 

L (m+1 ) 


o6 


ITl 




- hp m 

§_■) 

2L ^ 


( 2 . 19 ) 


For 



and 



2.0 

CL 


The ectuarion (2.10 redi. ces to the foilovang dimensionless 
f orra 



3 L J U 


m-1 

m 


( 2 . 20 ) 


The solution of the equation (2.20) is represented 

in figure (9), "by plotting < 7 ^ versus on semilogarithmiG 

% 

paper. From the plot it can be seen that discharge ratio 
increases with increasin.' value of “X and m. For all values 

U 

of greater than 1 . 0 , c^m greater than 1 , 0 , i.e. nonlinear 
discharge is mor^ than linear discharge. V/hen h) equals to 

U 

1 . 0 , the effect of la on discharge ratio is nil, 

2 , 3 . Flow into a trench - Combined Artesian and Gravity Flow 

2 . 3.1 Figure (3) illustrates a trench in a stratum under 
combined artesian and gravity flo^/» 

Let Ijq^ be the distance from the face of the trench 
to the point at which flow changes from artesian to gravity 



type . I’or value of y upto flo'" la gravity type and 
for y between and I flow is aitesxan type. 
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for aiteFian flo'' 


V 


JL 

BI) 


Uoinp enuation (1 .3) 


1 = 


dh 

dy 




in 


( 2 , 21 ) 


Integra cing equation (2.21) for boundary conditions at 
y = Ig, h = D 


h-D 

(y-Lo) 


n / Q \in. 


( 2 . 22 ) 


At y = I,, h = H 


n/'-Q \m 


/ H-B 


) = 


0 


(2.23) 


v/ith Q = 


when m = 1 , Q, = BDK . ( ?^) 

Iquation (2.^ can be reduced to the following form 


(2.24) 




rii 



m-1 


5 


- 1=0 


(2.25) 
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where ^ 




2 or f lo’.. , 


~T _ 


Q 


1 


_ dh 



m 

) 


( 2 . 26 ) 


Incecrating the equation (2.26) for h with boundary 
conditions at y = 0, h = hg and y = Lq, h = B. 


B“+1 - h”'^'’ 


in+1 


= 0 1, 


(2.27) 


v/ith 0 = ^ 


^ra+^ __ ^^m+1 



m+l 


m 

l^oC . 


,m-1 


- h^ m 

( —) . Lr 


(2.28) 


2 B, 


G 


The value of discharge (Q) is same in equations 
(2.25) and (2.27) by continuity, 

Bistance 1^ at which flow changes from gravity to 
ortesian type is found out from equations (2,23) and (2,27) « 



iQuatin^ the equstionr. (2,23) and (2.27) for the value 
ox Q , 


2G 


'vher 



L 


Por 


when 


X • 0 • 


H-D 

(L^ 


Iv 


(n+1 ) 


h 


2 r 


d2 


(2.29) 


D(1-r^2^) 


m+l 


(H-D) (m+1) + D d-rJJ') 


(2,30) 


H 

^ = 2.0, equation ( 230 ) simplifies to 


L 


( . m+1 V 


(ra +2 - ) 


(2.31) 


ra =1 , 


(D^ - h^) 


(2DH - D^ - hf) 


same as linear case. 


(2.32) 


The equation (2,31) is represented in figure (10) 

^G 

versus m. The ratio ^ 


with values of 


versus m. The ratio 


decreases with 
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iKcreasi.-, ■’'■alTjes of . 1 . The efXect x& laoie ax low values 
of 1,3. upbo 0.25. For value of tne effect 

of r Gi. ratio lo insignifiv-^'rt . 


2,4 Radial Plo\; Into a ’'el] Completely Penetrating a 
Confined Aquifer (Artesian Plow) 

2.4.1 Case (I) : 

The problem under consideration is one of radial 

sjc'-ixc xr^'" and in illustrated in figure (11), The flow is 

a 

assui’ed to be t\'o dimensional and S(.quifer is homogenous 
and i&otropic. The area of flow at a distance r from 
the center of the well is 2 77 rh, where J> is the thickness 
of aquifer. 


V = 


_Q 

2 iTfD 


1 


^h 

dr 


aQ 

2 7rrD 


+ 



4 TT 


(2.37) 


since at r = r^^^ , h = h^ 


h 


aQ 

h + 

2TrP 


loge (“^ — ) + 

V/ 



1 

r 


) 


ai r R, h = H 


(2.34) 



(2.35) 


aQ p 

H = h + loge (-^) + 

2 tT"D 'V 






Ihe discharge ratio iB obtameri froTi 

- 1 =0 

ioge 

vhere = r^j and ^ ) 

The drawdora curve can be expressed as 


( 2 . 36 ) 


h - hg ^loge (r/r,,) (l-od) (r-1) 

H - \ ’ lose (r^j) (r^ 3 -l ) "''^5 ^ r 


The discharge ratio calculated from equations (2.36) 

IS plotted as c ;6 versus 'V on somi logarithmic paper In 

6 

figure ( 13 ) j with ratio of radii of the third parameter! 

This discharge ratio (cy(j ) in bhis case also decreases 

logarithmically v/ith increasing value. The reduction 

in dj is more for the range of 'y) between 0.1 and 1,0 

^ 6 

than 1,0 to 10.0, Increasing r^^ value also have decreasing 
effect on dj , However for a tenfold increase in ^(^ 3 ? 
value of 06 decreases by about 0 , 1 , and r^^ above 100 , the 
decrease in (dj is further reduced. 
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2,4.'"'' Case (ii) - Power La,/ Colution 

In thifc ca&pj solution t .■ tho above problem is 
obtaiued by equation ( 1 .3) 

1 = cv”^ 

dh 

dr 

0 

V = 

2 TTrD 


At r = h = 


h-h^ = C (- 


Q 


m 


r 


2 rr D 


( 1 -m) r 


m~1 




At r = R, h = H, 


H - h =0 ( — 2 — )“ 

® 2 rE 


d-m) 


•%m-1 


J.T 


r. 


m-1 


w 


The discharge ratio cL = is ootaxned from 

m-1 ^ 


oC - - h 


) “ ( E“-^ - ) 5 lose (-^) 

V/ 




m-1 


( 1 - (~ — )“ loge r^j 

^d3 ) ^ 


(2.38) 


(2.39) 


m 


(2,40) 
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where ^ 


7 


CR 


and 



R 


1 


\v 


The solution to t?iO ’’ouftion - rc] resented in 
figure (H) for as function of ^7 . The incj'easing 

values of ^ and increases tno dischcxge ratio [ oC ) 

logcri'chmically and naxiruuin dischaige rate for a given value 
of ^ and r^j occurs v/hen m = 2.0, The discharge ratio 


incred&es considerably for 
for lower values. 


oi 

Lrj 


beb^feen 1,0 and 10.0 than 


2.5.1 Radial Flow Into a Well CoHpletely Penetrating 
an Unconfined Aquifer ( Gravity Flow) 


Case ( 1 ) 


An equation for steady radial flow to a well 
in an unconfined aquifer using Frrchheircer ' s equation (1,2). 
Referring to the figure (12^ 


dh _ Q.a _____ 

(ir 2rrrh 2 , 

4 rr h' 


(2,41) 


where 2 tj" rh is the area of flow at distance r from the 
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center ol the v'ell . 


2.2 dh 




I? 


aQ 

2rr 


rh - 


y2l 


= 0 


(2.42) 


Let 


r 


= R. 


h 


h. 


= H. 


dr 


he 

r„ iE^ 


(2.43) 


Substituting equation (2.43) in (2.42) 


dH. 


dR, 


aQ 


bQ‘ 


R.jxH.1 


2Tr^ 4TT^(rw'hJ) H^aR^ 


(2.44) 


If ^ *5 (dimensionless discharge) and 


% 

h 


= r 


d4* 


dHn 


dR, 


Q* 


Q 


2R.,^H.| 


4 2 


■d4 


R^x 


(2.45) 
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J'or linear cage equation can be written in the 
follo\/iiig form 




Qd^, . loge (-t: — ) + 1 .0 




(2.4G) 


and Qd^. is diraensionlesr' linear discharge and is equal 

t. 

1 0 


hllTK 


The solution to the equation (2.45) is represented 
in figures (15) (16) (17) and (18), Since the exact 
solution IS not possible for the equation (2,45), 
numerical solution vras obtained by Runga-Kutta method, with 
the aa d of digital computer facility. 


In the figure (15) the ratio is plotted 

r ® 

againsL — p- for different values of Q*, keeping parameter 

av 

and r^^ constant. The solution to fche equation (2.46) 

is also reuresented in this figure for comparing linear 

and non~linear flows, For the same discharge (dimensionless) 

and the effect of nonlinearity increases the 

dra\/dov/n ratio (|^ — ) . Por linear case increase in ^ 

e B 

value is gradual, Por exaniple for the discharge =0,1, 

for increase in ~ value from 1.0 to 10.0, drawdown 

ratio increases by about 10% for linear case whereas for 

nonlinear case, drav/down ratio increases by about 10 times 

i1:g initial value. 
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h 


In the figure (16), ^ versus 


— IS represented 
’e ~w 

for sare range of Qv, , '-nth hiH*lier values of^ i.e, 

^2 

nonlineerity parameter. In this case also the plot is 
similar bo previous case, Ho\/Qver at lo\/er value of • 

betv^oen 1.0 and 3.0, affect of nonlinearity parameter ) 

h ^ 

on T- 


ratio IS pronounced. 

Figure (17) shov/s the plot of 


h 


h. 


versus ~ for 
^v; 


values of'I^ ranging from 0.001 to 100.0 for some value 

of Q^<. . The effect of nonlinearity parameter (*7 ) on 

2 

drav/down ratio is insignificant for value of 0) upto 

e ^2 

10.0. For '7) equals to 100, drawdov/n ratio increases by 
'■2 

about ,05 . 


Figure (18) shows the plot of 


h. 


versus —zt' for 
w 


varying values of r^^, keeping 'T’ and constant. From 
the figure (18) it is evident that r^^ value has no effect 
on drawdown ratio, since for all values of r^^ drawdown 
rapio IS represented by single curve, 

2,5,2 Casa (11) ; Power law Solution 
From equation (1.3), 


dh 

■5r 


m 


CT 


= G (- 


Q 


m 


2 


•) 


(2,47) 


Integrating the equation (2.46) for h with boundary 
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condj-tions , the following equa'cion (2.43) can be obtained. 


“_1 ( 1 -r-j) ( 21T ll?. 


1 +in 


^ - c TT^ 


- h 


1 


(h'-" - r;-“) 


(2.48) 


'/hen m = 1, Q = Qj linear case. 


\ / 


ibh h = Qj^ and = 


R 


, the equation (2.48) may 


'v/ 


be reduced to The following form. 


o(: _-(§ ) 

^ {!«) L 


v/here 

^8 4i 


1 

Tj 1 +ra m ii-1 

(g-) -1. 


r''“^ - 1 

^d3 ^ 


1 


m 


g loge (2 


The solution of the equation (2.49) is shown in 

figure (19), for ck! as function of with ra as third 

® / 

parameter. For values of m upto 1,2, increase in c(^ for 

all values of ^ is gradual i.e, ^versus ^0 relationship 
‘‘6 

IS alno&t a straight line. For all values of r^^ and 
less than 0.10, variation in is insignificant and oC 


.49) 


increases with increase in the value of m. 
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FIG 19 UNCONFINED WELL FLOW (POWER LAW) 
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chapter pit 

n'L'UPP or rc’T-LiiEAR pi.o'Y 

3»1 Purpo&p of the Experiment 

Thiy invest igdtioa deals with the study of flow 
of v/cibor Liirough porous tnedia and epecificelly the 
ooef Lciei'bs ixi the Forchheiroerc eq.uation, 

3.2 Experimental Set-up 

A Schematic diagram of the experimental arrangement 
IS oho\/n 111 figure (20). A continuous flovv supply was 
provided by a overhead tanl: v'hicb maintainb constant 
head, The permeameter v/as 15 erns diameter tube and 
37.5 ems 111 length v/iih number of pieesure tappings 
alon^ the lengtn of the tuoe. To prevent tne channeling 
along the cylinder \,'a31s, the pressure taps were placed 
in tyliG permeameter so that oaoh Lap v/as 90° around the 
cyli icier Crom the tap belov'. A thin metal plate vath 
perf orotioiis was placed be] 0 '' the permeameter to support 
the material filled in the best section. The details 
of the permeameter tube is shov^n in the figure (21). 

The permeameter tube was connected bo a long pipe, 

15 ems in diameter, which in turn v/as connected to the 



5 '!- 


ovcrhoJu'' Lc'iik. viilvc j roviclca ot thp junction 
Goiuicc fc unj;, the vertical pi e section to the lon^ pipe 
at olic .groLind level, contiolo the flo'' rate through 
c]ip porMcamet er . The w'Hter alo< through tiie r'ermeameter 
verliccll/T- do\/n\^rardG , valve pi-ovided at the downstream 
e':d of the permearnet ere connect! che 15 erne pipe with a 

7.5 (Iierneter pi | e wliich rests inincu' the measuiin^ 

t. nh. The mes&uiing tank, 90 ers x 75 ems, '/hu-ch collects 
the watei ut the exit end ol the pipe \/as provided v/i th 
a 90'^ ^^-notch to measure the cliscliar^'^ a dr^in adjacent 
to the jj. ensuring tank collects the \ axer ' hioh in turn 
f]o\/ to a sump and a centrifugal i)ump lifts the water 
to cho overhead tank to complete tho cycle. Pressure 
head at various tappings were measured v/ith a mercury 
manoxtieter, Anothei manoxieter was used to check the 
Roasuromonts , 

5.5 Procedure 

ITaterials - In this investigation 3 types of 
mfitt I'lair were used, rai/iely river gravel, granite gravel 
and £lass spheres. The physical properties of these 
mabericJs are given in the table (5). 

VoT the studies of river gravel the permeameter 
v/as filled with the river gravel of median diameter (d50) 



55 


0«44 cmo, Ti^uro (22) "cc 3 hei^'ht of 50 gt.i& '''ithout 
the i|iateriai. fho density uao D*Gint9i2icd 
unifori by placing fcbc na:;eriol at coiistant height, 

A i.crecn (Ko. 54) placed on tue perlcrateci - etel 
plaLe and also at Lho tapping point c to prevent choking, 
'The pori'ieai icter was properly firceJ to the siipply pipe 
tind joirl's v/ere filled v/ith white lecd to prevent leakage, 
Vr]vo‘'’ at the upstream end Oiid at the do''n&treara end 
were opened and water v/as allowed bo flo'./ through the 
uedunn. The material v/as thoroughly saturated and air 
bul)b!'''c at che end of the tappings were carefully removed. 
The pj essure head at various tappings v/ere measured v/ith 
a mercury ipcnometer. 

dfter caiibruting the y-notch, discharge through 
V~notcli was measured with the aid of a point gauge. The 
bej-ivorature during experiment w'as noted. The procedure 
was repeated for various heaas and corresponding variables 
were 2 "ea&ured . For river gravel the experiment v/as 
iop('aled Tor 6 runs with different porosities. 

The grante gravel used in the experiment were 
of two sizes; oi median diameters (d 50) .547 ems and 
,5 13 cma with porosities .45 and ,465 respectively. 
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In th" run rnT'^bei 9, Lnher-s':. oi 

rii . cL.' 'fGi^e uc;eci c'-j n _ '-On’ra 'vns ii lied 

oI jC cnf’ ah teXcre a'id 'I'ea’-jielar v'-ig 
vit-’iced uairn c vibrator to I'e^p tiie poroe 3 ty 
ilia i e] ’1 a 3 uni f orir . 


1 .5 cms 
to a 

t •■oroughl^ 
oi fche 


A total of r runF \/itli different luatenals 
foa ■'i'3xj oue values of j-oj oeities 'vere conducted and the 
aiUiiiii-- r/ of i:he experiiient '1 daia are given in tables 
(5) -'Ud ( 6 ), 


5*4 rLCsulte and Diocussions 

The figure (23) represents the plot of ^ versus 
V foi different runs. There exists a stright line 
re] at Lonship for this plot and values of coefficients 
\/ere obtained from this plot- The table ( 5 ) shows the 
val.uos of a and h nlitr med for 9 runs for different materials, 

fig. (24) sho'-s the ilot of friction factor (f) 
versus foynolds number (f^) ^'iffe'rent runs. In the 
same figure two typical cur^'-es on - foi river gravel 
(J[^ ) and another for river &and (X^) '^l&o shown. 

The curves (X^ ) and (X^) are for Dudgeons data (Lee-1968), 
lu IS interesting to note that the experimental curves 
lies between the two curves (X^ ) and (X 2 ) for measured 
iteynolds number (R^) upto 1500. Por higher Reynolds 
nuril.or (R ) range measured the experimental curve is 
sliglicly shifted downv'ards. 
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'■Thus fro'u the relation ::hip between friction 
fecboc (P) and Reynoldj number (P J it can bo concluded 
for ie;;/no]ds number ran^e 10~200'j. Porehheirnera equation 

Lk- Vet j id, 

fpeciiic perTioabilitY, 1: =K— — 

\diere K= Darcy's coefficienc oL permrabiiity 
= Dynamic viscosity 
^ - Spocific weighs of fluid 

fho specific permeability (k) is used while 
coi L "Ic cing a and b values. 

The plot of specific permeability k versus 
copf'hcjcnt a ig shown in figure (25) for the reported - 
data. Tt is inlererting to note a very good correlation 
exinbs bebvifeen the values of k and 'a' and the slope 
oL the line is 45°. The relationship between k and a 
may be written in the following equation as 

ak = 10"^ (3.1) 

The values of 'a’ obtained in the present 
bbudies 'were plotted against k in Pig. (26). A straight 
line correlation on a log-log plot is seen. The range 
of Reynolds number 200 -1 500, However when the line (X^) 
of figure (25) was plotted on this figure (26), it is 
seen bhat the line (X^) is shifted to the right. This can 
be abbributod as the effect of Reynolds number (Bg). 
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'■Thus it can be concluded that a and k, and the 
propel bi onality constant is a fu’ otion of Reynolds 
nunbei The good correlation of Figure (25)? can 

be uOed 10 predict the coefficient a for known k if the 
Reynolds number is in the range 10-200, 

Figure (27) sho\/o the plot of specific per- 
meability Ic versus b for the reported data, In this 
ca.ge also a good oorrelaticn exists bet'/zeen the values 
of k and b. The relationship between k and b is given by 
bhe equation of the form 

bk = 1 .42 X 10“^ Cb 2,) 

The values of 'b* obtained in the present studies were 
plobted against k in Figure (28). A straight line 
corrolabion on a log-log plot is seen. The range of 
Reunion d& number 200-1500. Ko’^evor v/hen the line (X^) 
of figure (27) was plotted on this figure (28) , 

15 seen that tlie line (X^) is shift ea lo the right. 

This can be attributed as the effect of Reynolds number (H^). 

Thus it can be conoluded that b and k and 
proporGionality constant is a function of Reynolds 
number (Rg). The good correlation between k and b, 

Figure (27), can be used to predict the coefficient b 
for known k if the Reynolds number is in the range 10-200, 
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ProPictLon of Valuer, a and b ; 

■'or a l^no'"r iiearly imifonfi I'arLiclo size, "the 
f^necifnc permeability in a functiOxi of particle size 
and till' variation is shown in fi^^ure (29). Hence usmg 
using Figure (29) (Fur'er - 19C9), alcng \/ith the figures 
(?5) and (27), the coefficients a and b dh non-linear 
flo\, through a granular material of FeynoldxS number (Hg) 
range 10-200 can be predicted, 

jicaiaplc - For a particle size of 5 in,m. from the plot of 
spcciiic peimeability versus particle size, Figure (29), 
ipociiac permeability value is 2 x 10 ^ ems^ . fhe value 
of a and b can be obtained from the relationship of 
k versus a and k versus b, Figures (25) amd (27)» For 
pari, jcle rj 3 ze of 5 ni.m. and specific permeability 
2 X 10*"'^ cm. 3 ^ the vali>es of a and b are ,05 and ,041 
respectively. 




FrG-20 SCHEMATIC DIAGRAM OF THE EXPERIMENTAL SET-UP 
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FIG. 22 GRAIN SIZE DISTRIBUTION OF RIVER 
GRAVEL 
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Fr<5.24 VARIATION OF FRICTION FACTOR WITH REYNOLDS NUMBER 
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FIG 29 RELATION OF SPECIFIC PERMEABILITY 

TO PARTICLE SIZE FOR NEARLY UNIFORivl 
POROUS MEDIA (RUMER.R R 1969) 
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FIG. 30 'NON HOMOGENEOUS POROUS MEDIA . 
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chatter IV 


KOTR iIHLaA pro,/ THRCUGH HOIT-HOHOGEFOUS TCROjo 


FliDIA 


4 t 1 Hie flow oV v/£cter through porouo media is usually 
Obudiod as^unmg that the media is ivsotrpic and 
ho2iJO, jenono , Often this lo not; true where gravel and 
cooroe sard of different sizes are used in layers as 
in filter beds or grovel packs of tube wells. In such 
coses non-horaogerfous tiedium can be convened into an 
cquivolont single homogenous and isotropic layer for 
Lhe pui'pose of analysis. It is advantegeous if the 
equivalent nonlinear permeability coeflicients (a,b) 
oi bhc I’edinm are itrown. Similar relationships are 
avO'ilable for harcy's law v/herein equivalent coefficient's 
of permeability for la-''-erer'i soil can be easily obtained 
(Hari - 1962 ). 

Figure (30) lilustx'ater a vertical section 
through a stratified porous media ox H tiieir isotropic 
layers ol thickness d-| , d^ • ..... with nonlinear 
coefficients (a^ 9 b-| ) 9 ^2^ **’ ^^n^ ^n^ respectively* 

4.2 Analysis 

(a) Plow Parallel to Stratification 

For flow in the direction parallel to stratification 
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the clicchirge thr'ough each layer is governed by the following 
equetiono. 


1-1 = 


i_ ~ 




= a^V^+b-jV^ = 



= UpV^+b^V^ = 

^2 

Up 



0 

= «nVVn = 

®'n 

(-■) 


Q, 2 


Qo 2 

^2 


' 1 ' 

^2 


(4.1) 


Qv, 2 




n 


Alco Tor horizontal flow 


ll ^ lo 


= 1 


n 


\/'heio 1^, I2 ot;c- hydrarlic gradient in successive layers 

for an equivalent single homogenous layer 
of bui ah. esr < 1 ', the hydraulic grodienb ‘1' is 


1 - 1 1 - I2 - . . . • ij^ 


= 


, 0 , 


/W V /Q V 2 

= (^,) + (^,) 


where = 


+ d^ 


+ <3.. 


n 


(4.2) 

(4.3) 
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Addin: 1' eqnctions e'^UHx^on (^.1) leaoo to 




NO M 




(4.4) 


Oubstibuting Q^ = = ra^Q, .. ' ('^•5) 


+ iHg + . . . m^ =-' 1 .0 



1 .0 


(4.6) 


Subcbitutinp for Q2 .... equation (4.4) and 

equribin^' the coerfic3ente ot equal }io''rr of Q bo zero. 



”!e 

1 



(4.7) 


and 





(4,8) 


ii’or obtaimns the values of , lUg ... the 
cquabion (4,5) is substituted in equation (4,1) to give 
bhe follow3 2ig relationships. 



+ 


h| 

— 9 ni^2 

d, - 


’2 ^2 




a 


+ 


» 2 \ 
rA = A 

“ (4.9) 


4 * m 
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In equation (4.9), , nir, , '^r expressed 

jii terms of A and oilier known qi nbitie-^ ar 1 A is 
debcrniined from the relationship 


"n 


1 


Knowing values of rn^ , ... 

deLermined from equations 


the value F. of and 

(4.7) a 'd (4.8) respectively. 


(b) Flow Normal to Stratification * 

In this case, the discharge per ■' nit ar^a of 

each layer and total head loss throuih layers is equal 

to sum of the head losses in each layer . 

N 

li' = ^ 3n 
1 

\/hen 1 = a^y -h byV^ (4.11) 

and lj^= + bj^V^ (4.12) 

oubstituting the above relationships of i f^nd and 
equating the coefficients of equal pov/ers of Y to zero, 
gives following relationships, 

A ^ (% «n)/ ( £ , -In) 


(4.13) 



11 

V = 5 ,^ i-„ V d^) 
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(4.14) 


4.3 Experijvenbal ''Verification - f’lov/ tTornal to Stratifica bion 

Eyperiiiient y were condiicted on two layers of 
^ravol packed in the perineanet er . The size of coarser 
oravo] used was 5.47 m.m. { d^Q ) and finer gravel of 
size 3.168 m.m. (^^50^ thickness and the materials 
'k'Cx'p placed ab consiant height with densitiee^l ,375 gws/ 

CO and 1 ,305 gms/cc respectively. After satuiation, 
experiiiionbs \/erc conducted by varying the head over the 
per ii'00.,ie ter and measuring the gradienb and the di&cuarge. 

The suimrarj^ of the data are tabulated in the tables (7) a-nd 
(Ca and 8b) , 

4.4 rierAilla and Discussion : 

Tht head lose obtained Cor 4 different heads 
snou ohat there is a good agreement between the head- 
loss obtained by using the equations 4.10, 4.11 and 
measured values. The measured head loss and calculated 
hr>' d loss are nearly same, the percentage of error 
in al] the cases are less than 5%, v/hich sho\/s the 
accuracy with which the above equations can be used to 
estimate the head loss in layered porous media. 




Sii the samo e't’jatiori cc i b- u'-cJ to predict the 
hoad losr. bhrG..^h any nui/oej lay-- red >’oJia \/ith known 
valucG of GoefLicieiit fox’ each layer, 

4«5 HcOu loss Thrcu^.l Gra’/el fackr loi Tube 7ell& : 

The floe of water tbioeyl rr "el pacKs for tube 
Vi^ells is governed by non-linear flov' equations. Tt i& 
possible to estimate the head loss through grcivel pack 
(ti]Ler' pack) for known discharge and nsek thickness using 
[’03 ehheimers equate on, 

For a tube well drilled in eixhcr typo of 
aquifer (uniform or nonuniform) and shrouded by 
unifo3’iJi gravel designed on the basis of C.B.I,P. 

3''ecoiiiflenda bions (Garg - 1970) sand free discharge is 
larjUi'od fur gravel pack thickness of about 12.5 ems. 


4*6 Analysis for Plow Through Tv/o layers of Packing ; 

let the thtcknecs of packing be t^ and ^ 2 . 

T]k' rndjus of tube well is and 1-2 be the radius 

oC the coarser and fiiier gravel pack respectively. 




'■'m = + ( *1 + ^2^/2 


where mean raciiup oi uh: _.rr^^el pack 



— 5 I 

By conLinuiiy t^ischar.'^e through each layer in same 

■'' = 2tt = 2 tTr^2 (4.13) 

'/jic c ai'’ the T'ean velocities at mean radius 

and respectively. 

Head loss through grave] pack layers I and II 
Ua 0 g LVi n by 

■■•hi = *1 

ipdp = (apVjj^p + ^ 2^2 ^ ^2 


(4.16) 



Ascsuraing flo\. is rcidial and norpial to stratification, 


equivalent non-linear coefi icienL s lor t’/o lajers c?n be 
1 OLiiicl as follown. 

Tobal head loss id' 

= 

where a , are equivalent nonlinear coefficients and 
iy I 

V ’tea.' velocity of flov/ for entire section. Comparing 
m 

I, he cocLLicients ol same poi/er, 


v>« h h 

(aiY„,)t| t (agV^g) tg 
“y' (ti + tg) 

, _ *g (4.19) 

UiiiU I b - 

v/ (t, + tg) 


Since 


Q 


2 TT ^m 1 


V. 


Q 

2 TT r^2 


V. 


m 


Q 

2Trr^ 
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( ^ + ^ ) 


‘m I 




a = 

y 




( 4 . 20 ) 


axid by = ( 


^m1 


+ 


^:n2 


)/ 


(b^+ig 


)A 


ni 


(4.21 ) 


ir = -12 


Then a = 

t/ 


fm2 





( 4 . 22 ) 


(4.23) 


Knaving the values ol the equivalent 
hi.ireiOie gradient (i) and heed loss through entire 
.yuloM can be found from the relaxion 


ia'= (a^V^ + byV^) (t^+tj) 


(4.24) 



CHAPTER V 


7a 


5 1 Conolusions 

In this investigation, many field situations 
\/here in non-linear type of flow may occur are high 
lighLed, Solutions to some of the non-linear flow 
problems are presented using Porchheimers and Missbachs 
eq,uatioiis. The solution of these problems indicate that 
non-linear type of flow, discharge virill be less than the 
discharge under Darcy flow condition. For a given 
discharge drav/down will be more in case of nonlinear 
flow , 

The experimental investigation shoves that there 
1,3 a good correlation between the nonlinear coefficients 
and b and specific permeability Ic. This correlation 
can be used to predict the values of a and b for nearly 
unilorm particles with Reynolds number (Eg) range 10-200, 
The friction factor (F) versus Reynolds number (Rg) can 
be belter represented by different lines for various 
materials over wide range of Reynolds number (Eg). 

The equations are developed for obtaining 
equivalent nonlinear coefficients for flows both along 
and normal to the direction of stratification of an 
W-layered porous media. The experimental verification 
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Ghow.ci "tliGi; bhere is good, correlaiion bet\/een the measured 
head loss and the head loss obtained from the equations 

(4.10) ai id (4.11), 

Hov/ever further analyfcicol and experimental 
otudiCG ore necessary Lo the solution.^ of the problems 
involving higher fiov/ velocities and non homogenous 
po''''fJur. jiied-Sj such as flo’o tlirough rock fill dams and 
Li] Ler bedj. 
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APPEIv'DIX - A 


table 1 


(DATA OP DUDDEOII-1964) 


SI. 

a 

^ 2 / 2 

Medium 

Stand 

Specific 

Particle 

Wo. 

sec ./cmo 

secf/cm 


ard 

Peiueab- 

Size in 





error 

ilifcyp 

ems . 





V 

k Dm. 


1 

2 

3 

4 

5 

6 

7 

1 

0.0115 

0.0162 

Blue 

Mebal 

.0020 

1 .03x10"’^ 

1 .1 

2 

0.0073 

0 .0077 

Marble 

,0006 

1 .75x10'^ 

1 ,58 




Mixture 



3 

8.1161 

-0.0961 

Wepean 

Sand 

2,0478 

1 , 66x10"^ 

0.027 

4 

0.1904 

0.2174 

River 

Gravel 

.0029 

7.07x10“^ 

0.2 

5 

0.7891 

0.2232 

River 

Gravel 

.2873 

1 .69x10'^ 

0.095 

6 

0.1 661 

0.095 

-Blue 







Metal 

.1184 

8.05x10 ^ 

0.19 

7 

0.0779 

0 .0573 

Blue 

Metal 

.034 

1 .72x10''^ 

0.47 

8 

0.0024 

0 .0051 

River 

Gravel 

.0005 

5 .4x10"^ 

4.0 

9 

0.0035 

0.0121 

Blue 

Metal 

.0008 

4.0x10”^ 

2.50 

10 

0.0082 

0.0145 

River 


1 .62x10"^ 

1.2 



Gravel 

.0032 

11 

0.0189 

0.0262 

Blue 

Metal 

.0067 

7.0x10'“'^ 

1.9 

12 

0.0061 

0.01 17 

Blue 

Metal 

.0018 

2.2xl0~^ 

1 .82 


IX 


1 

2 

3 

4 5 

6 

7 

12 

0.0061 

0.0117 

Blue ,0018 

Ibtal 

2.2x10"^ 

1 .8 

I3 

0 .0143 

0.22 

Blue .0034 

Mebal 

9.2x10"^ 

1 .05 

U 

0.01 1 

0.0103 

Marolos .0017 

1 .2x10“^ 

1 .56 

15 

0.0058 

0.0066 

Marblec .0042 

2. 24x1 O”^ 

2.46 

16 

0.00064 

0.0015 

River ,0003 

Gravel 

2.47x10“^ 

1 .56x10”^ 

8.4 

17 

0.005 

0 .0063 

Marblec .0001 

1 .56 

18 

0.0076 

0.0095 

Marbles .00018 

1 .73xlO~^ 

1 .56 

19 

0.0051 

0 .0097 

Blue .0007 

Metal 

2.3x10“^ 

3.18x10’"^ 

1 .05 

20 

0 ,0036 

0.0049 

Marbles .0014 

2.85 



ill 


TABIE 2 

(LATA OF AffilEL fi CJIMA-I g69) 


ax. 
Fo . 

Env6'& Eigat; or 

a 

SCO ./cmo 

^ 2 ? 
aec . /cm 

Hed luri 

Gtsnd 

ord 

error 

S 

y 

Specific 

Permeab- 

ility 

2 

k cm. 

Par Cl cl e 
eizo in 
ems . 

1 

2 

3 

4 

5 

6 

7 

8 

1 

Ahmed 

7.59 

0.745 

Sand 

1 .5 

2.1x10“^ 

0.054 

2 

/ihmcd 

3.8 

0.454 

oond 

1 .7 

3.96x10"^ 

0.0764 

3 

Ahmed 

2.3 

0.308 

Sand 

1 .0 

6.91x10’’^ 

0.107 

4 

Alimed 

1 .49 

0.24 

Sand 

2.5 

1 ,0xl0~^ 

0.14 

5 

Alimecl 

0.938 

0.179 

Sand 

1 .0 

1 .69x10“^ 

0.199 

6 

Alimod 

0.694 

0.165 

Sand 

2,0 

2.21x10~^ 

0 ,258 

7 

Allen 

1 .47 

0.142 

Oi'anulsr 

3.3 

8 .6x10"^ 

0 ,0855 

8 

KLalro 

0.149 

0.0623 

Absorbent 

Glasc 

1 ,6 

6 ,7x10’'^ 

0 .32 

9 

Bro\mell 

0.0647 

0.0183 

Beads 

Glass 

1 .4 

1 .5x10'“'^ 

0.53 

10 

Brownell 

0.089 

0.021 

Beads 

Uickel 

5.0 

1 .1 2x1o”^ 

0.35^ 





Saddles 






iv 


3 4 5 6 7 8 


11 

Prancher 

16.6 

7.96 

Ottcwc 

S'ind 

4.1 

8.2x10 

0 .07 

12 

Porchhirn'^r 

0.00408 

0 .0005 

Sand 

2.2 

2.3x10“^ 

0.5 

13 

li'or ciihimer 

0.0123 

0 .00092 

Sand 

2.1 

7 .6xl0“''' 

0.3 

U 

Kirkhrm 

0.00895 

0.0117 

Rbrble 

4.1 

1 .19x10-^ 

1 .6 

15 

Tiiiidquist 

0.0674 

0.0368 

Sand 

3.0 

1 .38x10“^ 

0.492 

16 

Gj-nclnuiat 

1 .164 

0.292 

Sand 

4.5 

8.0x10*“^ 

0 .105 

1? 

Mc'liasheri 

0.189 

0.137 

Sand 

6.5 

4.94x10"^ 

0.5 

18 

oun 'do 

0 .1 45 

0.064 

Glass 

Spheres 

4.4 

6.45x10“^ 

0.3 


Y 


TABLE 3 


(BATA OF Ryvl-IGAITADHA BAO AITB &URESH-1970) 


SI. 

No. 

a 

sec ./cms 

b p p 

Sect/cnY 

Mecl lum 

Porosity 

(n) 

Speciixc 

Permeab- 

ility 

k cm.^' 

Particle 
size in 

cms . 

1 

2 

3 

4 

5 

6 

7 

1 

0.099 

0.263 

River 

Gravel 

Pailirly 

Round 

.40 

7.30x10”^ 

0.101 

2 

1 .15 

0.345 

rt 

.381 

6.53x10"'^ 

0.101 

3 

0.325 

0.11 

II 

.436 

22.54x10“^ 

0.17 

4 

0.475 

0.199 

11 

.417 

15.90x10"^ 

0.17 

5 

0.400 

0.164 

II 

.403 

18.78x10'"^ 

0.17 

6 

0.515 

0 .333 

II 

.392 

14.88x10"^ 

0.17 

7 

0.135 

0.072 

IT 

.430 

56 ,65x10*"^ 

0.286 

8 

0.225 

0.088 

11 

.423 

34.60x10“^ 

0.286 

9 

0.540 

0.400 

II 

.403 

22.10x10"^ 

0.286 



VD 


12 3 4 


10 

0.075 

0.053 

River 

Gr'^'vel 

Pairly 

Round 

11 

o 

o 

0.078 

t( 

12 

0.043 

0,043 

tl 

13 

0.074 

0.055 

n 

14 

0.105 

0.078 

IT 


5 

6 

7 

.304 

9.79x10"^ 

0.404 


.367 

68.47x10"^ 

0.404 

-372 

1G5 -8x10-^ 

0.55 

.456 

102.80x10’'*^ 

0.55 

.346 

73 .38x10“^ 

0.55 



Y3 1 


TABLE 4 

(BATA OF AMDAKMCHI'-r & VABDAE/uJLU-1 963 ) 


SI. 

Fo. 

Medium 

in 

1/G = (FJ) 

C .G.S .Unit p 

1 

2 

3 

4 

1 

Co^a-se Sand 

1 .11 

4.14x10'^ 

2 

rr 

1 .3 

2.43x10~^ 

3 

w 

1.63 

2.44x10"^ 

4 

Medium Sand 

1 .07 

3.12x10'^ 

5 

tl 

1 .26 

1 .89x10"^ 

6 

tl 

1 .48 

1 .75x10“^ 

7 

?ine Sand 

1 ,11 

0.55x10’’^ 

8 

w 

1 ,60 

0.069x10“^ 

9 

Very Pine Sand 

1 .00 

0.170x10"^ 

10 

t! 

1 .07 

0.147x10"^ 



VilX 


table 5 



BHYG 

I cat proj 

ExITIEG OP 

POPCDS 

I'lSDIA 



Bun 

No. 

Mefiium a 

secs./cms sec&./cric^ 

Porosity Specific 
Peri'iea b- 
ility 

2 

k cm. 

Particle Reynolcis 
size number 

(cI^q) Run^e 

ems ^e 

1 

2 

3 

4 

5 

6 

7 

8 

1 

River Gravel 

0.3076 

0.0855 

.40 

1 .72x10“^ 

0.44 

452-890 

2 

T1 

0.675 

0.10 

.59 

1 .785x10"’^ 

0.44 

275-807 

3 

n 

1 .096 

0.0961 

.363 

1 .67x10"'^ 

0.44 

212-627 

4 

n 

1 .267 

0.058 

.386 

1 .83x10“'^ 

0.44 

333-770 

5 

It 

1 .465 

0.067 

.394 

1 ,905x10”'^ 

0.44 

212-633 

6 

u 

0.800 

0.1125 

.366 

1 .72x:10'’'^ 

0.44 

317-663 

7 

Crranito 

CtUci vel 

0.57 

0.06 

.465 

2.24x10~*^ 

0.318 

272-732 

3 

tt 

0.283 

0.0916 

.45 

2.3x10”'^ 

0,57 

657-1549 

9 

Glofjo 

Spheres 

0.105 

0.007 

.396 


1 ,5 

4500-7577 



TABLE 6 


SUMMARY or THE EXPERII'IIirAL RESULTS 


SI. 

Ho. 

Run 

ITo. 

Q 

Lib/Sec . 

V 

Cms ,/Sec . 

1 

Ee 

(TO ) 

^ ^ .. ■ ■ ■ 

2) 

p 

Sgijo’- 

1 

2 

3 

4 

5 

6 

7 

1 

1 

1 .989 

10.54 

15.70 

653.75 

122.0 

2 

t! 

2.31 

12.24 

20.10 

759.2 

115.8 

3 

n 

2.462 

13.05 

22.40 

809.44 

113.55 

4 

IT 

2.58 

13.67 

25.20 

847.9 

116.42 

5 

It 

2.709 

14.35 

27.00 

990.07 

113.2 

6 

tl 

1 .505 

7.97 

10.2 

494.35 

138.62 

7 

11 

1 .81 

9.59 

13.50 

594.83 

126.72 

8 

II 

2.072 

10.97 

17.20 

620.42 

123.38 

9 

Tl 

1 ,376 

7.29 

8.48 

452.17 

137.75 

10 

It 

2.604 

13.80 

23.80 

855,96 

107.89 




1 2 3 4 5 6 


11 

2 

2.556 

12 

ti 

2.386 

13 

11 

2.20 

14 

tr 

2.003 

15 ' 

It 

1 ,849 

16 

If 

1 .66 

17 

M 

1 .534 

18 

n 

1 .322 

19 

11 

1 .067 

20 

11 

0.872 


13.54 

25.20 

12.64 

23 .00 

1 1 .66 

21 .40 

10.61 

18.60 

9,80 

16 ,50 

8,79 

13.80 

8.31 

12.50 

7.00 

9.80 

5.65 

7.24 

4-62 

5.16 


807.16 

118. 66 

753.50 

124.27 

695 .09 

135 .88 

632,49 

142.63 

584 . 20 

148,30 

524.00 

154.18 

484.70 

163.22 

417.30 

172.65 

336.81 

195.79 

275.41 

201 .70 





1 

2 

3 

4 

21 

3 

2.10 

11 .12 

22 

ri 

1 ,915 

10.148 

23 

rr 

1 .778 

9.42 

24 

f[ 

1 .593 

8.44 

25 

n 

1 .517 

8.03 

26 

11 

1 .239 

6 ,56 

27 

rr 

1 .426 

10.55 

28 

If 

1 .139 

6.03 

29 

ri 

0.712 

3.77 

30 

n 

0.980 

5.19 


5 6 7 


23.80 

627.90 

166 .16 

21 .50 

573.02 

180.23 

19. 2 D 

531 .91 

186 .79 

17.30 

476.57 

209.66 

15.50 

453.42 

206.1? 

12,30 

370.42 

246 .74 

13.90 

426.32 

210.51 

9.97 

340.49 

236.71 

4.67 

212.87 

283 .65 

CD 

« 

O 

293 .06 

259.59 



XI 1 


1 

n 

3 

4 

5 

6 

7 

31 

4 

2.03 

10.76 

20.30 

625 .49 

151 .36 

32 

4 

1 .843 

9.76 

17.90 

567.36 

162.22 

33 

tl 

1 .722 

9.12 

16 .00 

530.16 

166.06 

34 

tt 

1 .523 

8.07 

13.75 

469.12 

182.26 

35 

tl 

1.393 

7.38 

12.70 

429.00 

201 .30 

36 

11 

1 .081 

5.73 

9.20 

.-Kl 

OJ 

O 

• 

o 

241 .89 

37 

IT 

1 .213 

6.43 

10.80 

373.78 

225.50 

38 

II 

2.128 

11.28 

21 .10 

655,72 

143.15 

39 

M 

2.229 

11 .81 

25.20 

686.53 

143.59 

40 


2.370 

12.56 

26.30 

730.13 

143.92 



XJ-13_ 


1 

2 

3 

41 

5 

1 .355 

42 

11 

2.143 

43 

TT 

2.051 

44 

f / 

1 .895 

45 

H 

1 .778 

46 

\i 

1 .183 

47 

It 

0.998 

48 

11 

0.719 

49 

tl 

1 .691 

50 

rt 

1 .376 


4 5 


7.18 

13.00 

11 ,35 

25.20 

10.87 

22.80 

10.04 

21 .40 

9.42 

20.60 

6,27 

11 .90 

5.29 

9.30 

3.81 

5.07 

8.96 

18.90 

7.29 

13.80 


6 

7 

400-48 

217.69 

633.07 

168.87 

606.30 

166.58 

560.00 

183.27 

525.42 

200.41 

349.72 

261 .51 

295 .06 

286 .89 

212.51 

501 .51 

499.76 

203.23 

406.61 

224.17 


XIV 


1 

2 

3 

4 

5 

6 

7 

51 

6 

1 .053 

5.579 

7.95 

317.78 

220 .49 

52 ^ 

It 

1 .529 

8.10 

13.60 

461 .32 

178.94 

53 

tr 

1 ,66 

8.79 

15.45 

506.79 

172.62 

54 

n 

1 .778 

9.42 

17.40 

536 .56 

169.27 

55 

ir 

2.093 

11 .09 

22.40 

631 .68 

157.23 

56 

II 

2.20 

11 .657 

25.40 

663.98 

161 .56 

57 

fi 

1 .942 

10.29 

20.70 

586.11 

168.76 

58 

tr 

1 ,817 

9.627 

19.10 

548.35 

177.91 

59 

n 

1 .239 

6.56 

10.00 

373.65 

200.60 

60 

If 

1 .076 

5.70 

5.70 

524,67 

151 .45 



XV 


3 4 5 6 7 


61 

7 

2.10 

62 

li 

3.395 

63 

t1 

3.166 

64 

If 

2.962 

65 

11 

2.784 

66 

II 

2.588 

67 

II 

2.348 

68 

f! 

1 .836 

69 

11 

i .599 

70 

ft 

1 .264 


11 .13 

10.90 

17.99 

28,20 

16.79 

26.30 

15.69 

24.10 

14.75 

21 ,40 

13.71 

19.45 

12.43 

16,50 

9.73 

11 .05 

8.475 

9.45 

6.70 

6,44 


453 .00 

54.20 

732.21 

53.68 

683.36 

57.47 

638.59 

60.31 

600.33 

60.59 

558.00 

63.74 

505.91 

60.79 

396 ,02 

71 .90 

344.90 

81 .05 

272.69 

88.38 


xvx 


1 

2 

3 

4 

5 

6 

7 

71 

8 

2.923 

15.51 

18.40 

1038.0 

82.08 

72 

ir 

2.229 

11.81 

10.25 

866.57 

78.86 

73 

ft 

1 .691 

8.96 

4.81 

657.45 

64.30 

74 

IE 

2.517 

13.34 

15.70 

978 .84 

94.68 

75 

TT 

3.442 

18.24 

26.00 

1338.0 

83.87 

76 

It 

3.148 

16.68 

23.30 

1223.0 

89.87 

77 

t? 

3.704 

19.62 

29.80 

1439.0 

83 .08 

78 

U 

3.987 

21 .12 

34,20 

1549.0 

82.28 

79 

tt 

3.834 

20.31 

31 ,50 

1 490 .0 

81 .95 

80 

It 

3.509 

18.59 

26.30 

1364.0 

81 ,67 



XV3 1 


1 

2 

5 

4 

5 

6 

7 

81 

9 

7 .111 

57.67 

13.50 

7216 

28 .00 

82 

T1 

7.467 

39.56 

14.60 

7578 

27.45 

83 

11 

8.006 

42.41 

17.30 

8123 

28.30 

84 

n 

7.741 

41 .01 

15.15 

7855 

26,51 

85 

It 

6.85 

36.29 

12.80 

6951 

28.60 

86 

ti 

6.47 

34.28 

11 .90 

6566 

29.80 

87 

ti 

5.971 

31 .63 

11 .05 

6058 

32.50 

88 

II 

5.724 

30.32 

9.45 

5808 

30.25 

89 

11 

5.371 

28.45 

8.62 

5450 

31 .34 

90 

11 

4.436 

o 

C\J 

6,44 

4501 

34.32 



XVill 


TABLE 7 


iiiY.rtCfa ncTi^vi'Li cr tfi; ■'v' 


IiLiteriBi Layer I PoroGity(n) Layer II Porosity 


GrciiiLo Size 


Size 


Gi^Tve] 4.7-6 25 
mm 


1 .656-4.7 
mm 


.505 




R 

e 


table a (b) 


iJUi I A];Y C" Tjji' K 


1' T 


va 


/ ^ 2 , ^ ^ 2 , 

sec./ secv sec./ secT/ 

cms cms'^ cius. cits. 


1 < A. 

Calcu 

Meas Perce 

la tod 

ured ntage 

Head 

Head error 

loss 

Loss 

in cii.s. 

in cms , 


1564 

21 .1 

0.869 

0,0143 

1 .025 

0.0145 

425 

432 

1 .75 

1338 

20,6 

0.869 

0.0143 

1 .025 

0.0143 

415 

399 

3,86 

1178 

17.25 

0.765 

0.0201 

0,965 

0,0178 

354 

339 

4.25 

1254 

19,00 

0,765 

0.0201 

0.965 

0.0178 

382 

389 

1 .95 


1 .95 


APPEL]I;IX B 


HEAD LOSS THI^OUGH GRAVEL PACK OP A TUBE '/ffiLL 
In thiG appendix, the c£ Iculntion of head loss 
fchrou;’h ''grovel pack of a tube ’uell using the equation 
(4.2'1-) iL illustrated with a sir^ple example. 

Bor a tube well of diameter 15 cms. v/ith 
gra'^'e! pack thickness of 12.5 cms and aquifer particle 
size of 0.5 m.m. (medium sond)^ the discharge is 20 
litres per metre length of the screen Pig. (32). 

In this case, equation (4.24) simplifies to 
Id' = + bv2 


v/here V is the mean velocity at radial distance r. 


m 



'C|r/^ 


Figure “32 


r^ = Mean radius of the gravel pack. 

Discharge = Average area of flow x mean velocity. 


20x10' 


= 2 TP ^ 27 , 5 X 100 V 


m 


V =2.31 cmfl'/sec. 
m 

Assuming pack aquifer ratio as 12.0 (Garg - 1970). 
Gravel pack size = 12x0. 5 = 6 m.m. 

Value of coefficients a and b for gravel pack of 6 m.m. 
size are assumed as .065 &nd .042 respectively. 
Substituting the values of a, b and in equation, 

Head loss = 4.52 cms. 



xi:i. 


I :: c 

, I ..i.'icT: li' :?*! ■ T 'I r. 

2he Lahora-uox-y i, £,hoils oi the i >6 nre abilities 
iae-- L f iiienx r> of pnrcuf i r. cari’ief* ovt in 

pCi fc" - L 0 e r 01 txnitc cx’Ob' acctioii. 'I'he ic 

bo7’-i" ■' y 2 oiiooch \/all ’/hicn la i. pGj j G'ble. The 

I’x't.ei^co c_ aif'ooth 'vall ^fXocte lL^ x'Lo’' behi.i''’-iour 

fco the »/all './Inch in lunri oxiGcta the velocity 
xoi a thiii zone near the '.'03 1 i& highei then that for inner 
zoxic nninl,; due to the highei roroeicy near xhe wall 
z on -■> . 

Tn Llie prerent Bnaly&i^ The f] o / is sosumed to 
Lowin’, linear ]a\; in the inner zone and for \fall zone 
vjJocity IS hirlif^r the flo'w is governed by non- 
n -Cir ] Dvdgeon (1967) eatmated the wall effect 

i: rri ...oter w'lth. so'oe afesumpticnn and using the 
power lo" ecuaticn (l,3)« 

figure (33) shows cross section ot a permeamet er , 



Figure 33 





Ili. jI' I L* Ci-C > iTi'ictei (.1 permecii eter and t, vid 
' -J' ::rr.e 

, - O-lTi r_rc.i''iiol re-, of inncx’ Z'^xife (D-2t) 

d_:‘i ^ter 

“ . ot al .'i.r ea 




A 


\/ 1 

= Zone Area 


Fcx’ In lor Zone 


K 


.( 01 ) 


V 

L 

ro-- 


1 


= Ve] ocxty ixi tlio inner zone 
’.'all zone, I’emf, e^'ua ti on ( I .'i') 


a V + b 

M V/ V' 17 


(C2) 


' iiero c,.. b, are noniineai coexfiojente lor jail 
, ve. ocic;," ihiouj.h vain zone. Oince hydraulic 
LO Jjxneo and ooter zones ore same. 


? 


JL 

Z 


= a V J b, 1 / 

'7 W '7 V/ 


b j- a V - 

\/ V/ 'J \7 


Lh oi 

in 


zone and 
gradient 

(C.4) 


C 


(C.5) 



XZClli. 


colvixi Liic qusdr^tic equation (C.5) for end taking 

pc"'ioi 'C 

(G.6) 

2t)2 

(C.7) 

Ciibsti Luting equation (0.6) in (C-7) 
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